Electron-hole asymmetry in magnetic properties of lightly doped high-Tc 
superconductors: a quantum Monte Carlo study 
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Using a recently developed variational quantum Monte Carlo method, magnetic properties of 
high-Tc superconductors are studied at zero temperature (T), by numerical simulations on the 2D 
t-J model. Our focus here is to explore the difference in the properties of p-type and n-type cuprates 
as a function of the carrier concentrations close to half filling. As observed experimentally, it is found 
that the antiferromagnetically ordered phase persists even for a small, yet finite amount of carrier 
doping, and that this phase is more robust in the electron doped case. 



There have been stimulating studies on high-Tc super- 
conductors to explore the underlying mechanism and un- 
usual normal state properties. Experimentally much at- 
tention has been recently devoted to rather lightly doped 
p- and 71- type cuprates [l| . This is partially because high 
quality samples have become now available, and because 
it is now clear that understanding the evolution of the 
electronic state with carrier doping from undoped Mott 
insulator is crucial to clarify the nature of the systems. 
Theoretically it is still in the middle of extensive debate. 
It is widely accepted however that a 2D i- J- like model 
can describe low-energy physics of cuprates defined by 
the following Hamiltonian 2] ; 



H = (S. • S, - n.n,/4) - t ^ (^^S,- . + H.c, 



(1) 



Here the notation is conventional (see e.g., Ref In- 
troduction of next nearest-neighbor hopping t' deserves 
some explanation; small cluster calculations have shown 
that the t-J model indeed capture low energy physics 
of more involved multi-band models with t' /t negative 
(positive) for hole-doped (electron-doped) case, at least 
for small carrier densities x 0. We remind that this 
sign difference originates from the fact that the parent 
cuprates are charge transfer insulators, and therefore a 
"carrier" represents in the model a Zhang-Rice singlet 
(an inert Cu configuration) for the hole (electron) 
doped system. It is also now widely accepted that, in 
the half-filling undoped case, the model shows antiferro- 
magnetic (AFM) long-range order at T = 0. It is very 
interesting therefore to study theoretically how strong 
this AFM phase is against carrier doping and how the 
robustness of the phase depends on the nature of doped 
carriers. 

There exist several numerical studies on small clusters 
(< 26 sites) indicating that at finite x close to half-filling 
positive (negative) t' /t enhances (degrades) AFM corre- 
lations compared to t' = 0]. Although these results 
might give the correct trend at short distances, it is dif- 
ficult within such small cluster calculations to really as- 



sess the existence of long-range order. For this purpose 
we have studied this issue using quantum Monte Carlo 
techniques, which allows much larger sizes (~ 100 sites). 
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FIG. 1: Absolute value of real-space spin correlation functions 
for 2D t-t'-J model with J/t = 0.3 and t' /t = -0.2 (top), 0.0 
(middle), and 0.2 (bottom). Carrier densities x are 0.04 (left) 
and 0.08 (right) on a tilted square lattice of \/98 x \/98 . Open 
squares and solid circles are for VMC and FN, respectively. 



To study ground sate properties of the 2D t-t'-J model 
Q , we employ a recently developed variational quantum 
Monte Carlo method by Sorella 0] . Results presented 
here are computed by variational and fixed node methods 
(VMC and FN, respectively), and the details on how to 
minimize the energy for a variational wave function with 
many parameters are found in Ref. j^. We just mention 
here that we employ the most general singlet variational 
wave function of the Jastrow-BCS form, with a pairing 
function of d-wave symmetry. The Jastrow factor, which 
includes the Gutzwiller projector, makes the wave func- 
tion of the RVB type. After all the possible variational 



2 



parameters are optimized, we obtain a fairly good guid- 
ing function that is used for FN. Also note that the vari- 
ational energy of FN is lower than the one of VMC . 
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FIG. 2; Same as in Fig. [3 but a; = 0.12 (left) and 0.16 (right). 

Main results are summarized in Figs. ^ and [3 where 
real-space spin-spin correlation functions are given for 
various t'/t and doped carrier concentrations x 0| at 
J/t — 0.3. These results provide clear evidence that 
strong spin correlations remain finite at the largest dis- 
tances for finite carrier concentrations x. The critical 
Xc above which long range order disappears apparently 
depends on the parameters used; for the present set of 
t'/t = -0.2, 0, and 0.2, Xc ~ 0.06, 0.08, and 0.12, respec- 
tively. Of particular interest is the result that the spin 



correlations are enhanced (degraded) when positive (neg- 
ative) t'/t is incorporated into the model, correspond- 
ing to electron (hole) doing. Our results are consistent 
with the previous studies 0, 1^ and experimental obser- 
vations 

In conclusion, magnetic properties of the lightly doped 
2D t-J model with next nearest neighbor hopping t' are 
numerically studied to understand the different behaviors 
as a function of carrier concentrations for hole {t'/t < 0) 
and electron {t'/t > 0) doping. Our large cluster cal- 
culations clearly show that AFM long-range order exists 
away from half-filling and is more robust for the electron 
doped case. These results compare very well with exper- 
imental observations. More detailed analyses and other 
properties will be presented in a separate paper [i(tI |. 
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